Collective modes of a two-dimensional spin-1/2 Fermi gas in a harmonic trap 
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We derive analytical expressions for the frequency and damping of the lowest collective modes of 
a two-dimensional Fermi gas using kinetic theory. For strong coupling, we furthermore show that 
pairing correlations overcompensate the effects of Pauli blocking on the collision rate for a large 
range of temperatures, resulting in a rate which is larger than that of a classical gas. Our results 
agree well with experimental data, and they recover the observed cross-over from collisionless to 
hydrodynamic behaviour with increasing coupling for the quadruple mode. Finally, we show that a 
trap anisotropy within the experimental bounds results in a damping of the breathing mode which 
is comparable to what is observed, even for a scale invariant system. 



I. INTRODUCTION 



A new generation of experiments realising two- 
dimensional (2D) atomic Fermi gases provide a unique 
possibility to systematically explore 2D quantum sys- 
tems, including phenomena relevant to high temperature 
and organic superconductors, and 3 He films [TH3] ■ The 
study of collective modes has proven to be a powerful 
way to probe the properties of 3D atomic gases includ- 
ing the effects of pairing and strong correlations [SHE|. 
The damping of these modes is related to the transport 
properties of the gases, and in particular the damping 
of the quadrupole mode is controlled by the shear vis- 
cosity in the hydrodynamic limit. The shear viscosity of 
quantum Fermi fluids has received lots of interest due to 
the conjectured minimum of the ratio of viscosity and 
entropy [5141 lj. Recently, the first experiments concern- 
ing the collective modes of a 2D Fermi gas have been 
performed [12] . The breathing mode frequency was re- 
ported to be very close to that of an ideal gas with a low 
and constant damping throughout the different interac- 
tion regimes, which indicates a scale invariant system 
with a very small bulk viscosity [13H16) . The quadrupole 
mode was shown to exhibit a clear transition between 
collisionless and hydrodynamic behaviour with increas- 
ing coupling strength. Following this experiment, the 
shear viscosity for a 2D gas has been calculated using 
kinetic theory, and the result was then used to calcu- 
late the quadrupole mode damping in the hydrodynamic 
limit [TT1 - H5] . From the good agreement between experi- 
mental data for the quadrupole mode and a numerical 
solution of the Boltzmann equation, it was concluded 
that kinetic theory provides an accurate description in 
the temperature and coupling strength regime relevant 
to the experiment [20]. Here, we calculate the frequency 
and damping of the energetically lowest collective modes 
of a 2D Fermi g function of temperature and cou- 



pling strength from approximate analytic solutions to the 
Boltzmann equation obtained by expanding it in appro- 
priate basis functions. The effects of Pauli blocking and 
pairing correlations on the collision rate are systemat- 
ically analyzed, and we show that for strong coupling 
there is a large temperature regime, where pairing cor- 
relations dominate the effects of Pauli blocking. This 
indicates a large so-called pairing pseudo-gap region in 
the Bose-Einstein condensate (BEC) to Bardccn-Cooper- 
Schrieffer (BCS) phase diagram |2"lTf2l)] . Our predictions 
for the mode frequencies and damping are shown to agree 
well with the experimental results. In particular, the 
theory reproduces quite well the collisionless to hydrody- 
namic cross-over of the quadrupole mode with increasing 
coupling strength. By including the effects of a possible 
slight trap anisotropy, we show that the resulting cou- 
pling between quadrupole and breathings modes can lead 
to additional damping of the breathing mode even when 
the bulk viscosity is zero. This result can be important 
for a proper interpretation of experiments. 
Our paper is organized as follows: First we derive the col- 
lective mode frequencies using the method of moments. 
Then we proceed by calculating the collision rate using 
various levels of approximations. We demonstrate that 
the collision rate can be larger than that of a classical gas 
due to pairing correlations. Our results are then shown 
to provide a good quantitative description of the exper- 
imentally observed collective modes. Finally, we explore 
how a trap anisotropy within the experimental bounds 
results in a significant damping of the breathing mode. 

II. FORMALISM 

We consider a gas of N fermionic atoms of mass m 
trapped in a harmonic potential 

V(x, y, z) = -m (u 2 x x 2 + uj 2 y y 2 + lu 2 z z 2 ) . (1) 
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FIG. 1: (Color online) (a) Plot of damping T versus frequency ujqm for the quadrupole mode. The theory curve (solid line) only 
relies on the moment method and is independent of a microscopic model for the collision time r. The dots are the experimental 
data, (b) The collision rate 1/r as a function of interaction strength at T/Tf = 0.47. The inset shows the relaxation rate as 
a function of T at fixed interaction strength \n(kpa2Y>) = 0.5. All theory curves were calculated with (black) /without (blue) 
Pauli- blocking and with Pauli-blocking and medium corrections (red). The collision rate including medium corrections diverges 
when T* as given by the Thouless criterion (indicated by the vertical dashed line) is approached from above. 



There is an equal number N/2 of atoms in two internal 
spin states which wc denote j" and J., and atoms in differ- 
ent spin states interact via a short-range potential char- 
acterised by a scattering length a whereas the interaction 
between equal spin atoms can be neglected. We focus on 
the quasi-2D limit with u> x ,U! y <C u) z aim T,Ep <C lj z , 
where T is the temperature and Tp = Ep = Nui x L) y 
the Fermi temperature (For brevity, we set H = ks = 1)- 
In this limit, the motion of the atoms is frozen along the 
z-axis and the kinematics of the system is essentially 2D. 
We use the semi-classical Boltzmann equation 

(5 t +r-V r +p-V p ) / = -/[/] (2) 

to describe the collective motion of the gas, where /(r, p) 
is the distribution function, r = p/m and p = —V r V. 
Since we consider density modes, we have taken the same 
distribution function for the f and the J. atoms. As in 
Refs. [TTHT^] , we assume that the main effects of interac- 
tions are captured by the collision integral / so that inter- 
actions effects on the left side of the Boltzmann equation 
can be ignored [TH]. The collective modes can be ob- 
tained by linearizing the Boltzmann equation around the 
equilibrium distribution function 

/ (r,p) = e/3[ep+v(r) _ M] - 1 , (3) 

where j3 — l/T and e p = p 2 /2m. In order 
to linearize Eq. pL we write /(r, p) = /o( r , p) + 
/o(r, p) [1 - /o(r, pj] 0(r, p), where the factor f (1 - f ) 
is introduced for later convenience |27j . The linearized 
Boltzmann equation then takes the form 

/o(l - /o) (d t + r • V r + p ■ V p ) - -I[<p] (4) 



with the collision integral 

I[4>i] = I d 2 p 2 ^ J de'\T\ 2 {4>i + - h - 

/i/a(l-/s)(l-/4) (5) 

where <pi stands for <fi(r, pi) etc., and we have defined 
cPp = d 2 p / (2-k) 2 . The scattering of particles with mo- 
menta pi and P2 to momenta P3 and P4 with P4 = 
Pi + P2 P3 is described by the scattering matrix T. 
The relative mass is m r = m/2 and 6' is the angle be- 
tween the outgoing relative momentum p' r — (p 3 — p 4 )/2 
and the center of mass momentum P = pi + P2 . Note 
that P is conserved in the collision and that |p r | = \p' r \ 
with p r = (pi — P2)/2. In Eq. ^ and the following, all 
distribution functions are of the equilibrium form given 
by Eq. (§. 

In vacuum, a collision between a f and a J, atom with 
momentum p r and kinetic energy e = p 2 /2m r in the 
center of mass frame is described by the low energy 2D 
T-matrix [2"g] 

%{e) = ~ . \ (6) 

v ' m r ln(e*/e) + in w 

where e* = f^/ma 2 ^ with the 2D scattering length 

a 2 D = Izy/n/B e ~^~l 2l ^ a % with l z = \j^Jmuj z and 
B = 0.905 [IH]- 

The collective modes are solutions to Eq. Q of the 
form (f>(r, p, t) = e~ ZUJt (f>(r, p) where Re(w) is the mode 
frequency and Im(w) is the damping rate. We expand 
4>(r, p) on a set of basis functions which form a closed 
set under application of the left hand side of Eq. Q, i.e. 
we write ^(r, p) = J2i c i4'i ( r , P)- Several previous studies 
have applied this method (or variations of it) successfully 
to harmonically trapped gases, see e.g. [2"7l I5DH55] . and 
here we apply it to a 2D Fermi gas. 
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FIG. 2: (Color online) (a)-(b): Quadrupole mode frequency (a) and damping (b) in units of ft± = (uj^Uy) 1 ^ 2 compared to the 
experimental data of Ref. [12] . In the gray shaded area, the normal state is unstable towards pairing (within BCS mean-field 
theory). We have added an offset of 0.1f2j_ to all damping theory curves to account for other sources of damping present in 
the experiment Ref. |12] , The insets show the eigenmode frequencies and damping rates for an anisotropic trap (uu y — 1.05 u> x ), 
which leads to a better agreement between theory and experimental data. 



III. ISOTROPIC TRAP 

The trap used in the experiment in Ref. :12 is nearly 
isotropic in the xy-plane, and we therefore first consider 
the case lo x = io y = 

For the breathing mode, we use the basis set 

0i = n 2 ± (x 2 + y 2 ), <f> 2 = tt±(xp x + yp y ), 3 = pI+pI- 

(7) 

Since this set is closed under application of the left hand 
side of Eq. Q and the collision integral vanishes by con- 
servation of energy and momentum, one finds an un- 
damped breathing mode at = 2il± [30l [34]. This 
is a consequence of harmonic trapping and rotation sym- 
metry, and the result agrees surprisingly well with what 
is seen in experiments |12j . 

For the quadrupole mode, we use the basis set 

0i = n 2 ± (x 2 -y 2 ), 02 =0,±(xp x -yp y ), 03 =p 2 x -p 2 . 

The corresponding equation determining the mode fre- 
quency is [BJ 



regimes. The effective collision rate 1/r which comes 
from the non-zero collision integral of 03, is 



1 



Jd 2 rd 2 p (p 2 x - p 2 y )I{p 2 x - p 2 y 



Jd 2 rd 2 pf (l-fo)(p 2 x 



(12) 



Before using Eq. (12) to calculate r as a func- 
tion of temperature and interaction strength, we 
first eliminate r from Eq. (|9| to obtain Tq = 

Im(w) is 



n 2 ± 



ujq where Tq 



the damping of the quadrupole mode with frequency 
ujq = Rc(w) [BJ. In Figfl] (a), we compare this with the 
experimental data of Ref. [T2] by plotting the damping 
rate as a function of the frequency. The good agreement 
indicates that our approach for calculating the collective 
modes captures the most important physics of the experi- 
ment. In particular, it shows that extending the basis set 
given in Eq. (|8| to include higher order functions [33j will 
most likely only yield a fairly modest gain in accuracy. 



IV. COLLISION RATE 



2 2 i / 2 2 \ n 



(9) 



where u> c \ = 2ft± is the collisionless and Whd = v2fii the 
hydrodynamic limit of the quadrupole mode frequency. 
Indeed, Eq. ^ has the solution 

U = Wcj - — (10) 

4r 

in the collisionless limit lot 3> 1, and the solution 

w = Whd - it. (11) 

in the hydrodynamic limit lot <C 1. This explicitly 
demonstrates the different roles of collisions in the two 



To make further connection to experiments, we must 
calculate r as a function of temperature and interaction 
strength. It is easiest to calculate r in the high tempera- 
ture classical limit z = exp(/3/x) <C 1, where we can ignore 
the Fermi-blocking factor (1 — fa)(l — /a) in Eq. ^ and 
use the Vacuum T-matrix given by Eq. Q. Equation 
( 12 1 then becomes 



1 

Tel 



TTZT 



G(e*/T), 



with the dimensionless integral 



G(y) 



dt 



o ln^(y/t)+^ 2 



(13) 



(14) 
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Converting the fugacity z to particle number using the 
ideal gas formula N = Nf + N± = 2T 2 z/uj x ujy gives 



1 



■kNlo x (jJ v 
AT 



G{e*/T). 



(15) 



Comparing this expression for the collision rate in a har- 
monically trapped gas in the classical limit with that for 
the shear viscosity ?/ c i of a uniform gas in the classical 
limit PUCES], we get 



T 



(16) 



where n(0) = mNoj x ujy/2irT is the density in the center 
of the trap. This is the usual relation between the col- 
lision time and the viscosity in the classical limit, with 
the factor 1/2 reflecting the trap average. From Eqs. ( 11 ) 
and ( 16 1, we furthermore see that the damping is propor- 



tional to the shear viscosity in the hydrodynamic limit 
as expected [i"7l \18\ [35] . 

For lower temperatures, the Pauli-blocking factor (1 — 
/s)(l— f±) in Eq. ^ for the collision integral significantly 
decreases the collision rate. On the other hand, medium 
effects in the scattering matrix were found to essentially 
cancel this reduction due to the pairing instability for a 
3D gas in the unitarity limit, resulting in a collision rate 
close to the classical prediction for temperatures all the 
way down to the critical temperature T c for superfluid- 
ity [BJ. In 2D, medium effects on the scattering matrix 
were even found to over-compensate the effects of Pauli 
blocking when calculating the shear viscosity and spin 
diffusion coefficient in the strong coupling limit [19] . In 
addition to keeping the Pauli-blocking terms in Eq. ^ 
for lower temperatures, we therefore also include medium 
effects in the T-matrix. This gives in the ladder approx- 
imation 136 



1 



1 



T M (P,e) T v (e) 



d 2 k- 



/k+P/2 + /k-P/2 

+ iO — k 2 /m — P 2 /Am 



(17) 



where iO denotes an infinitesimal positive imaginary part. 
Because center of mass and relative coordinates are cou- 
pled in the presence of a Fermi sea, this adds an extra de- 
pendence of the T-matrix on the center of mass momen- 
tum P. Fortunately, the numerical evaluation 1/r can 
be simplified significantly by using the symmetry/anti- 
symmetry of the collision operator under the exchange 
of its momenta, and the invariance under a simultaneous 
rotation of all momenta. Using this, the numerator of 
Eq. ( 12 ) can be simplified to 



4^2 



d 2 r / dPP / d 2 p r p 4 r / d9'\T MB \ 
Jo J Jo 



[1 - cos(2# - 2#')]/i/ 2 (l - /s)(l - U) (18) 

where 9 is the angle between p r and P. The denominator 
in ( 12 ) can be written as 



d 2 rd 2 ph{l~h)(p 2 x ~plf = - 



4T 4 m 2 



where Li s (a;) is the polylogarithm. The resulting expres- 
sion for the collision rate 1/r is evaluated numerically. 
We note that it only depends on T/T F because rescaling 
the x, y coordinates by the harmonic oscillator lengths 
l x _ y = 1/ ' -s/fnuJ^y cancels any additional dependence on 
trap frequencies. The fugacity z is found from to T/Tp 
via the formula for the non-interacting Fermi gas in a 
harmonic trap Tp/T = yj— 2Li2(— z). 

In Fig. [I] (b), we plot the collision rate as a function 
of interaction strength \n(kpa2r>) for T/Tp = 0.47 and 
as a function of T/Tp for \n(kpa2D) = 0.5 using three 
different approximations: The blue curves are the classi- 
cal expression given by Eq. |l3| , the black curves include 
Pauli blocking in the collision integral, and the red curves 
include Pauli blocking as well as the medium corrections 
to the T-matrix given by Eq. ( 17). The classical collision 



Li 3 (-z) (19) 



rate monotonically increases as the temperature is low- 
ered as a consequence of the increasing density. When 
Pauli-blocking is included this behavior is drastically dif- 
ferent, since the exclusion principle freezes out collisions 
in the degenerate regime T <^Tp. This freezing out re- 
sults in a maximum of the collision rate as a function of 
temperature as shown in the inset of Fig. [I] (b) . 

This picture is further modified when medium correc- 
tions in the T-matrix are included. Here the proxim- 
ity to the BCS instability towards pair formation gives 
rise to resonant behavior for particles colliding with op- 
posite momenta near the Fermi surface. In the strong 
coupling regime | \n(kpa2D)\ % 2, this causes a dramatic 
increase in the scattering rate over a large temperature 
range, diverging at the mean-field BCS transition tem- 
perature T* given by the Thouless criterion, as indicated 
by the dashed vertical lines in Fig. [I] (b). While there 
is strictly speaking no phase transition at T* in 2D, this 
should be considered as the temperature where pairing 
sets in, but fluctuations prevent superfluid order as long 
as T is above the Berezinskii-Kosterlitz- Thouless temper- 
ature Tbkt [261137] . The medium corrections in fact over- 
compensate the Pauli-blocking effects for strong coupling 
resulting in a collision rate which is larger than the classi- 
cal prediction, as can be seen from Fig. [l](b). This over- 
compensation is consistent with what was found for the 
shear viscosity and spin diffusion coefficient [19] . Medium 
corrections to the scattering in the strong coupling regime 
are thus even larger in 2D than in 3D [5J, which can be 
interpreted as the presence of a large so-called pairing 
pseudo-gap regime above the superfluid transition tem- 
perature Tbkt 37 . 

On the molecular side, beyond the BCS transition, the 
collision rate is again finite: In this limit, our theory 
effectively describes a repulsive Fermi gas in the upper 
branch [38] . Away from the BCS transition, the collision 
rate with medium effects approaches the Pauli-blocked 
result. 
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V. COMPARISON WITH EXPERIMENT 

We now compare our theoretical results with the ex- 
perimental data of Ref. [12] ■ It should be emphasized 
that the theoretical expressions for the frequency and 
damping contain no free parameters to fit theory and 
experiment. In Fig. [2] (a)-(b) we plot the frequency 
and damping of the quadrupole mode. As for the col- 
lision rate, we show three different curves: A classical 
calculation (blue), a calculation including Pauli blocking 
(black) , and a calculation where both Fermi blocking and 
medium effects on the T-matrix are included (red). We 
see that there is good agreement with the experimental 
results for quadrupole frequency and damping rates. In 
particular, the onset of the hydrodynamic regime around 
ln(fc^ a2_o) ^2 — 5 and the location of the maximum 
damping rate seen in the experiment agrees with our cal- 
culations. 

We see from Fig. [T] that the collision rate including 
medium corrections in the scattering matrix is in be- 
tween the classical result and the Pauli blocking result 
for ln(fc F a 2D ) > 1.6 and T/T F = 0.47. This explains 
why the frequency in Fig. [2] is in between the classical 
and the Pauli blocking prediction in this regime. When 
—0.6 < ln(fcp-a2D) 1-6 the medium collision rate is 
larger than the classical prediction, but in this regime 
the mode frequency is already very close to the hydrody- 
namic value, so there is no visible difference between the 
classical and the medium prediction for the frequency. 
The higher collision rate when medium corrections are 
included does however yield a smaller damping in the hy- 
drodynamic regime, as can be see from Fig. [2](b). Some- 
what surprisingly, it appears that including medium cor- 
rections gives a slightly worse agreement with the experi- 
mental results compared to the theory that includes only 
Pauli-blocking. Experimentally, when tuning from at- 
tractive interactions beyond \n(kpa2D) 1 one observes 
rapid heating of the gas, qualitatively in agreement with 
the expectation of increased pairing correlations yielding 
a higher three-body loss rate when entering the pairing 
pseudo-gap region T < T* . We note that our theoretical 
approach is not valid in the paired regime and further 
theoretical studies are needed to fully address the collec- 
tive modes at low temperatures. 



VI. ANISOTROPIC TRAP 

It is typically quite difficult to achieve perfect symmet- 
ric traps in experiments. In the Cambridge experiment, 
there is some uncertainty in the value of the trapping fre- 
quencies u x and uiy so that the trap may not be perfectly 
isotropic in the xy-plane. Instead, one typically has an 
upper bound of \u> x /ui v — 1| < 5% for the anisotropy. 
An anisotropic trap results in coupled breathing and 
quadrupole modes [30l 00] , and we shall now show that 
this can have important consequences on the interpreta- 
tion of the experimental results. 



We consider an anisotropic trap with lo x =^ lo v . To 
describe the coupled breathing and quadrupole modes, 
we use a generalized version of the basis functions Eqs. 

(ttT 



01 = u 2 x x 2 + u 2 y 2 , 4> 2 =uJ x xp x +uj y yp y , (20) 

2 



p x + P y , 4>4 =u x x -w y y , 

UJ x Xp x - UJyVPy, (j> 6 = P x - Pi- 



rn) 

(22) 



The calculation of the collective modes using this basis 
set is described in the Appendix. After a straightforward 
but rather lengthy calculation, we obtain that the normal 
mode frequencies are the solutions of 



(^ - w cll)(w 2 - Wcl2) + — (u 2 - ^hdlX^ 2 - w£ d2 ) = 

(23) 

with uj c \i/2 = 2^ x /y the coUisionless frequencies. The 
hydrodynamic frequencies are 



w hd,l/2 




14 



C0 2 x UJ 2 y +U y . (24) 



The collision time r appearing Eq. ( 23 ) is identical to 



the collision-time for an isotropic trap, when evaluated 



at the same fugacity 



since we can eliminate 



any dependence on trap frequencies by rescaling the real- 
space coordinates x,y. 

In Fig. [3] we show the frequencies and damping rates 
of these modes as a function of the collision time t 



for 



= 1.05 uj y . The coupling of the breathing and 



quadrupole modes is clearly visible. In the insets of 
Fig. [3j we compare our theory with the experimental 
data for the frequency and damping of breathing mode 
in Ref. [12]. The inset of Fig. 3 (a) shows the separately 
measured cloud width oscillation frequencies in x-, and 
y-direction of the breathing mode as a function of interac- 
tion strength. In the coUisionless regime (ln(fc F a2D) 1) 
the oscillations along the principal axis of the system are 
decoupled. The frequency difference directly reveals the 
system's anisotropy of about 5%. Towards the hydrody- 
namic regime the motion along the main axis start to 
couple leading to a single collective mode in the strong 
coupling limit (ln(fci?et2D) 3). The breathing mode fre- 
quency in that regime lies below the theoretical predic- 
tion, probably due to the anisotropy, which couples the 
monopolc mode with the quadrupole mode. An impor- 
tant consequence of this coupling is that there is a viscous 
damping of the breathing mode, which is undamped for 
an isotropic harmonic trap as the kinetic theory yields 
zero bulk viscosity. This potentially complicates the in- 
terpretation of the observed damping of the breathing 
mode in the Cambridge experiment, since it is necessary 
to separate the damping due to non-zero shear viscosity 
caused by a possible trap anisotropy from damping due 
to non-zero bulk viscosity. In fact, the inset of Fig. [3] (b) 
demonstrates that a small trap anisotropy of 5-10% can 
qualitatively explain the slight increase in damping of the 



1 Q ± t 2 3 4 1 n ± x 2 3 4 

FIG. 3: (Color online) Frequency (a) and damping (b) of the lowest collective modes of a two component Fermi gas in a slightly 
anisotropic quasi-2D harmonic trap (oj y = 1.Q5lo x ) normalized by the geometric mean of the frequencies S7x = {i^x^y) 1 ^ 2 
as a function of collision time r. The dashed lines in (a) denote the isotropic limits of the breathing and quadrupole mode 
frequencies ujb,q in the collisionless ojb,q/£}± = 2 and the hydrodynamic regimes ujb/£1± = 2,cjq/S1j_ = v2, respectively. The 
coupled breathing and quadrupole modes are shown in red and black respectively, and the blue curves correspond to the scissors 
mode (see appendix). Insets: (a) Decoupled eigenmodes in the collisionless regime (ln(fcFa2D) S> 10) start to lock towards the 
strongly-interacting regime (ln(fcFi2D) < 5), resulting in a single collective mode with strong breathing mode character, (b) 
Same plots for the parameters of the experiment |12| as a function of interaction strength at T = 0.47 Tf (here shown for 
the collision time obtained from the theory with Pauli-blocking) . The dots are the experimental data of the breathing mode 
damping rate for two different excitation strengths. [39] 



breathing mode |39j with coupling strength observed in 
Ref. [12] . This additional source of damping makes the 
bulk viscosity extracted from the damping even smaller, 
and more measurements are needed in order to pin down 
how small the bulk viscosity is and to compare with the- 
ory [HHE]. 
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VII. CONCLUSIONS 



We studied the collective modes of a 2D Fermi gas us- 
ing kinetic theory. Expanding the Boltzmann equation 
on basis functions, we obtained analytical results for the 
frequency and damping of the collective modes as a func- 
tion of temperature and coupling strength, which were 
shown to agree well with experimental data. We demon- 
strated that there is a large temperature range for strong 
coupling, where pairing correlations dominate the effects 
of Pauli blocking on the collision rate, resulting in a rate 
higher than the classical value. Finally, we showed that 
the coupling between the breathing and the quadrupole 
pole modes due to a slight anisotropy within experimen- 
tal bounds, may result in a damping of the breathing 
mode even when the bulk viscosity is zero. This can have 
important consequences for the interpretation of the ex- 
perimental results for the breathing mode. 
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Appendix A: Derivation of the equations for the collective modes in an anisotropic trap 



To find the collective modes, we expand the linearized Boltzmann-equation Eq. Q in a finite set of basis functions. 
This will allow us to obtain simple matrix equations for the collective mode frequencies and damping rates. Formally 
we are interested in solutions of the equation C u (j) = with the linear operator £ u defined as 



Cu>4> = /o(l - fo) + r • V r + p ■ V p ) 4> + I[<t>\. 



To solve this equation we project Eq. (Al) onto a set of functions </>(r, p) = Y^i a<j>i 
Boltzmann equation into a tractable finite dimensional matrix equation of the form 



(Al) 

(r,p), converting the linearized 

(A2) 



with coefficients given by 



J d 2 rd 2 p faCutfij 



f<Pr<Ppf Q (l-f Q )<ft- 
The collision integral I [(f)] vanishes for any function 

4>{r, p) = o(r) + p • b(r) + p 2 c(r) 



(A3) 



(A4) 



since it is local in space momentum and energy is conserved during collisions |27j . Applying this formalism to the 
moments Eqs. ( 20p2 ) yields then the set of equations for a general harmonic trap 



/ 


— iu> 


J (Wjc + UJy) 


ILJ 

2 







~(uj x + U)y) 




U) x + UJy 


UJy — UJ X 




iui 
2 


- \ (iO x + Uly) 


—iui 










2 ( W I _ Uy) 





— iui 




uJy ~ ^J x 





uJ x ~ t-^y 


~(UJ X + UJy) 


V 





2 fay — ^x ) 









UJ V 







I ("x 
I (^V ~ ^x ) 

I (UJ X + UJy) 

—icu 

- \ (y>x + u> v ) 





u> x — ui 



LJ X + Ul 

—iu) 



\ 



ii 



V 









C2 




C3 




C4 


V c 6 / 



(A5) 



This matrix equation has nontrivial solutions when the coefficient determinant vanishes. This results in Eq. ( 23 ) for 
the no rmal mode frequencies of the coupled monpole and quadrupole modes described by the basis functions Eqs. 
( 20|22 ). Here the only basis function for which the collision integral does not vanish is <j)e = p x — p y . By parity 
symmetry only the matrix element Mqq receives a contribution 1/r given by Eq. ( 12 ), giving a finite damping rate to 
the collective mode frequencies. For an isotropic trap, Eq. (A5) decouples into an independent undamped breathing 
at twice the trapping frequency [30] and a damped quadrupole mode with frequency and damping given by Eq. (j9). 
Finally, we mention that independent from the coupled quadrupole/breathing mode, there is also the scissors mode 
which can be described by the basis functions |41j 



f> 7 = uj x uj y xy, 8 = ui x p y x + uj y p x y 7 

i>9 = U x p y X - UJ y p x y, 010 = PxPy 



(A6) 
(A7) 



In the isotropic case ui x = u) y , this mode is simply one of the degenerate rotations of the quadrupole mode. The 
frequencies of the scissors mode are the solutions of the equation [4"Tl W2\ 



iuj 



(A8) 



where uj'. 



- 

mode 



hd 
U) x 



UJt 



uj 2 is the frequency of the scissors mode in the hydrodynamic limit, whereas ui' cll 



to,, and 



u)y\ are the mode frequencies in the collisionless limit. One obtains the matrix equation for the scissors 



— IbJ 

-\ (u) x +U)y) 

| fax ~ UJy) 



■ w, 







W X -r "Jy ^y 

—iui I { u x + uj y ) 

-iuj \{w x -u)y) 

-(0J x +U)y) U)y-UJ x ~IUJ + ^ J 




= 



(A9) 



8 



when projecting onto the basis functions Eq. (A6). Again this equation has nontrivial solutions when the determinant 



vanishes, leading to Eq. (A8|. In the isotropic limit, only three basis functions are coupled and in this case they 



simply describe the quadrupole mode [rotated by 45 degrees with respect to the quadrupole mode from Eq. (A5)] 
The collision rate for the scissors mode is given by 



J_ _ / d 2 rd 2 p p x pyl[p x py} 
t> ~ J d 2 rd 2 pf (l - k)plp 2 v 

Since both p x p y and p 2 — p 2 belong to the I — 2 representation of the rotation group, we have r' 
also valid in the more general anisotropic case where lj x ^ uj y . 



(A10) 



r. This result is 
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